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We compare the standard 50%-efRcient single beam splitter method for Bell-state measurement 
to a proposed 75%-eflicient auxiliary-photon-enhanced scheme [W. P. Grice, Phys. Rev. A 84, 

042331 (2011)] in light of realistic conditions. The two schemes are compared with consideration 
for high input state photon loss, auxiliary state photon loss, detector inefficiency and coupling 
loss, detector dark counts, and non-number-resolving detectors. We also analyze the two schemes 
when multiplexed arrays of non-number-resolving detectors are used. Furthermore, we explore 
the possibility of utilizing spontaneous parametric down-conversion as the auxiliary photon pair 
source required by the enhanced scheme. In these different cases, we determine the bounds on the 
detector parameters at which the enhanced scheme becomes superior to the standard scheme and 
describe the impact of the different imperfections on measurement success rate and discrimination 
fidelity. This is done using a combination of numeric and analytic techniques. For many of the 
cases discussed, the size of the Hilbert space and the number of measurement outcomes can be very 
large, which makes direct numerical solutions computationally costly. To alleviate this problem, 
all of our numerical computations are performed using pure states. This requires tracking the loss 
modes until measurement and treating dark counts as variations on measurement outcomes rather 
than modifications to the state itself. In addition, we provide approximate analytic expressions that 
illustrate the effect of different imperfections on the Bell-state analyzer quality. 


I. INTRODUCTION 

Bell-state measurements are a key component of quan¬ 
tum information processing protocols such as quantum 
teleportation [T], entanglement swapping [5] and dense 
coding They are particularly critical in the con¬ 
text of quantum communication. The implementation 
of large-scale quantum networks for long-distance quan¬ 
tum cryptography [J] and a future quantum internet 
requires the distribution of photon entanglement over 
long distances. Unfortunately, direct photon transmis¬ 
sion through hbers suffers from exponential loss over dis¬ 
tance causing rates to decrease below acceptable levels 
01 El [7]. In classical communication this problem is 
circumvented by signal amplihcation. However in the 
quantum world, this cannot be accomplished due to the 
no-cloning theorem for quantum states [SlEl. 

Many promising methods that overcome the direct 
transmission loss-rate trade-off use quantum repeaters 
[IMIl- These methods rely on generating entanglement 
over multiple shorter distances, and then connecting the 
links by entanglement swapping [2], to extend the en¬ 
tanglement over the entire distance. Free-space based 
proposals using satellites take advantage of reduced pho¬ 
ton absorption outside the atmosphere, allowing further 
distances [T8H22]- However, global entanglement distri¬ 
bution using satellites would also require entanglement 
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swapping due to earth curvature [22] . An ideal distribu¬ 
tion network may comprise a combination of fiber and 
free-space solutions, using fibers for regional distribution 
and satellites for intercontinental distribution. 

A shared requirement among all practical entangle¬ 
ment distribution protocols is successful entanglement 
swapping between neighboring nodes to create links. 
This is done by applying a projective Bell-state mea¬ 
surement (BSM) on two photons, one from each node 
0]. The measurement result would then indicate the 
entangled state of the link. Ideally, each measurement 
outcome should unambiguously indicate projection onto 
one of the four Bell states. Unfortunately, a perfect Bell- 
state analyzer cannot be constructed using only linear 
optical elements 02102] and cannot exceed 50% success 
rate without the addition of auxiliary photons |25j . The 
simplest Bell-state analyzer can resolve at most two of 
the four states resulting in a maximum 50% success rate 
for entanglement swapping |26| . This is a source of inef¬ 
ficiency in all repeater schemes and compounds with an 
increase in the number of swapping events. 

A perfectly complete Bell-state analyzer cannot be re¬ 
alized by simple means; however, efficiency can be in¬ 
creased arbitrarily close to unity at the cost of increased 
apparatus complexity and added auxiliary photons 07). 
As a consequence of increasing complexity, the design 
can become more susceptible to component flaws that re¬ 
duce efficiency and fidelity, offsetting initial gains. Hence 
it is necessary to compare BSM schemes in light of im¬ 
perfections such as non-number-resolving detectors, de¬ 
tector inefficiencies, dark counts, and photon loss rates. 
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To approach the problem of quantifying the analyzer 
quality, we use a combination of numeric and analytic 
techniques. Determining the exact quality of a Bell- 
state analyzer, when accounting for different imperfec¬ 
tions or auxiliary sources, can be analytically difficult 
and computationally costly. This is because the size of 
the Hilbert space and the number of measurement out¬ 
comes increases quickly when additional parameters, de¬ 
tectors, and photons are considered. Our approach to 
this problem is to use only pure states when perform¬ 
ing numerical computations, which reduces the required 
computational resources. This is accomplished by track¬ 
ing all loss modes until the measurement is performed 
and by treating dark counts as variations on measure¬ 
ment outcomes rather than altering the output state. We 
also perform analytic approximations to derive expres¬ 
sions that give insight into the influence of imperfections 
on the Bell-state analyzer quality and use our numerical 
results to verify the accuracy of these expressions. 

This paper is organized as follows. Section in¬ 
troduces the Bell measurement schemes, detector type 
models, loss considerations, and the alternative auxil¬ 
iary photon source used in our analysis. Section [m] out¬ 
lines the methods used to analyze the schemes under 
realistic conditions, including simple examples using the 
standard scheme as a case study. Numeric and analytic 
results are reported in section m which is divided into 
three parts. The first part explores replacing detectors 
with imperfect ones, the second considers multiplexed 
arrays of non-number-resolving detectors, and the third 
deals with an alternative auxiliary source. Finally, we 
summarize our conclusions in section 


II. BACKGROUND AND DEFINITIONS 
A. Bell-state measurement schemes 


A Bell-state analyzer is a device that performs a (par¬ 
tial) Bell-state measurement. It is characterized by its 
success rate for projecting an entangled state onto one of 
the four Bell states and identifying the result, assuming 
a uniform distribution of all four Bell states as input. 
These Bell states can be written as 


IV’'^)i2 = ^ (44 ± 44) |0)i2 
l'i^’'^)i2 = ^ (44 ±44) |o)i2 > 


where and w are photon creation operators for qubit 
mode and the subscript denotes spatial mode. 

The most standard 50%-efficient scheme for conduct¬ 
ing a BSM is composed of a single non-polarizing beam 
splitter and two detectors [26] , also referred to as a Hong- 
Ou-Mandel interferometer, as illustrated in FIG.(left). 
This is described by the beam splitter transformation 
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FIG. 1. (left) A schematic of the standard method for per¬ 
forming a Bell-state measurement, (right) Grice’s proposed 
first-order extension for Bell-state measurement. The spa¬ 
tial modes 1 and 2 represent the input pair while the spatial 
modes 3 and 4 represent the auxiliary pair. Gray bars rep¬ 
resent non-polarizing beam splitters. Diagrams are shown 
assuming measurement of time-bin encoded qubits. For po¬ 
larization encoded qubits, each detector in the above dia¬ 
grams is replaced with a polarizing beam-splitter and two 
detectors. For dual-rail encoded qubits, the above diagrams 
are duplicated, one for each rail. 


where is replaced by or U as appropriate. This Bell- 
state analyzer has at most a 50% success rate since only 
and of the four Bell states are distinguishable 
by measurement outcomes. The other two Bell states 
cannot be distinguished from each other but can be 
distinguished from The simple design and lack of 

auxiliary photons makes this method attractive for prac¬ 
tical implementations. However, limited success rate is 
undesirable especially when many measurements are re¬ 
quired to generate entanglement over long distances. 

We compare this standard scheme Bell-state analyzer 
to an extension proposed by W. P. Grice m- In his 
paper, Grice defined a series of enhanced Bell-state an¬ 
alyzers using increasingly complicated entangled auxil¬ 
iary states that, in the series limit, show an arbitrar¬ 
ily complete Bell measurement. The auxiliary states re¬ 
quired to approach a perfect Bell-state analyzer are im¬ 
practical to construct; however, the first-order extension 
only requires an auxiliary |^i >“'')34 Bell state to achieve a 
75%-efRcient Bell measurement, which is the enhanced 
scheme studied in this paper. This first-order extension 
utilizes four beam splitters and four detectors, as illus¬ 
trated in FIG. (right), and is described by the trans¬ 
formation 
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where spatial modes 1 and 2 contain the input state to be 
measured and spatial modes 3 and 4 contain the required 
auxiliary state. 
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B. Detector types 

Driven by the rapid expansion of quantum information 
science, single-photon detectors (SPDs) have seen signif¬ 
icant improvements over the past few decades [211 ESj. In 
particular, photon-number-resolving detectors (PNRDs) 
that use transition-edge sensors (TESs) have reached ex¬ 
cellent efficiencies of 95% at telecommunication wave¬ 
lengths [SOllH]. Unfortunately TES-based detectors op¬ 
erate at temperatures of 0.1 K and also have recovery 
time on the order of microseconds [22] ■ Long detector 
recovery time or ‘dead time’ is particularly undesired 
for applications that use time-bin encoded qubits, which 
are popular in quantum communication usi mm- 
Faster detectors such as single-photon avalanche diode 
(SPAD) detectors and recently developed superconduct¬ 
ing nanowire SPDs (SNSPDs) have recovery times on 
the order of nanoseconds [mnanz], which are more 
suited for use with time-bin qubits. Although they op¬ 
erate at higher temperatures (77 K - 250 K), SPAD de¬ 
tectors are relatively inefficient for telecommunication 
wavelengths (~ 10%) and cannot natively resolve photon 
number |2H1[22]- On the other hand, SNSPDs show im¬ 
pressive efficiencies of 93% for telecommunication wave¬ 
lengths jSH] but, like SPADs, they cannot resolve pho¬ 
ton number and, like TESs, require significant cooling 
(< 3 K) [22]. Currently there does not exist a detector 
type that satisfies all desired criteria for quantum com¬ 
munication applications; therefore careful analysis of the 
trade-offs is required on a per-application basis. In this 
paper we discuss and analyze some of these trade-offs in 
the context of a Bell-state analyzer. 

We consider four different types of abstract detec¬ 
tor models, two of which are relevant to time-bin en¬ 
coded applications only. The first and most ideal detec¬ 
tor model is a PNRD. In addition to resolving photon 
number, we assume a PNRD has negligible dead time. 
Secondly we consider a non-number-resolving SPD that 
cannot indicate photon number but again has negligi¬ 
ble dead time. That is, it is a binary detector (BD) that 
clicks when photons are incident and does not click when 
photons are not incident, similar to a SPAD detector or 
an SNSPD [221122] • We assume for simplicity that both 
the PNRD and BD can perfectly resolve qubit modes. 
This can be implemented, for example, by placing a po¬ 
larizing beam splitter prior to the detectors (for polar¬ 
ization encoded qubits), by placing a detector on each 
spatial rail (for dual-rail encoded qubits), or recording 
two time-separated events (for time-bin encoded qubits). 

The last two detectors we considered are similar to the 
PNRD and BD described above but have the additional 
flaw of long dead time. This flaw is only relevant to 
time-bin-encoded qubits provided the temporal modes 
are restricted to a maximum temporal separation. We 
describe a slow PNRD that can indicate photon num¬ 
ber but has a dead time longer than the time separation 
between qubit modes, which could occur with a TES de¬ 
tector [321ISI]. As a result, if one or more early photons 


are detected, the detector is unable to count possible late 
photons in the same spatial mode. A slow BD would be 
similar to a slow PNRD but again cannot resolve photon 
number, similar to a gated SPAD [2H1 [22] • To give an 
example, a slow PNRD measuring a state with two early 
photons and two late photons would indicate the pres¬ 
ence of two early photons whereas a slow BD measuring 
the same state would indicate only one early photon. 
In addition to number-resolving imperfections and long 
recovery times, each detector model is assumed to also 
suffer from detection inefficiencies and dark counts. 

Alongside the development of PNRDs such as the 
TES, detectors that attempt to resolve photon num¬ 
ber by multiplexing many non-PNRDs in either space 
[321EO] or time mi are also being pursued. These meth¬ 
ods attempt to combine the benefits of non-PNRDs such 
as warmer operation temperature or short dead times 
with number resolving capabilities. These multiplexed 
detectors can be modeled by arrays of imperfect detec¬ 
tors; however, undesirable effects caused by inefficiency 
and dark counts can be amplified by large detector ar¬ 
rays. We discuss regimes in which detector arrays of 
BDs and slow BDs can increase Bell-state analyzer suc¬ 
cess rates closer to those obtainable with PNRDs. 


C. Photon loss channels 

The input state that is to be measured can suffer from 
photon loss. In particular, when utilizing the Bell-state 
analyzer in repeater elementary links, the input pho¬ 
tons must first travel large distances. As a consequence, 
the input state to the Bell-state analyzer can have poor 
transmission rates m- Aside from a reduction in total 
efficiency of the measurement, this input loss can affect 
the fidelity of the scheme when considered along with 
detector dark counts. 

For the enhanced scheme, the auxiliary state can also 
suffer from photon loss. This loss channel represents 
a contribution from propagation absorption, coupling 
losses, and pair generation inefficiency. For bright aux¬ 
iliary pair sources, this loss rate is expected to be small 
as compared to the input state photon loss rate. How¬ 
ever, auxiliary state photon loss only affects the quality 
of the enhanced scheme, and not the standard scheme 
since it does not require auxiliary photons, and so this 
parameter is expected to directly affect the bounds on 
the regime where the enhanced scheme is superior to the 
standard scheme. 

Finally, both schemes are also subject to photon loss 
between the beam splitters and the detectors. This loss 
rate can be accounted for as combined with detector ef¬ 
ficiency to form an effective scheme loss/inefficiency pa¬ 
rameter. 
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D. Alternative auxiliary source 

Schemes that improve success rates beyond 50% re¬ 
quire auxiliary photon states. These states are suscep¬ 
tible to realistic imperfections. In some cases, using 
approximations for these auxiliary states may also be 
necessary for implementation to remain practical. The 
enhanced scheme studied in this paper requires an aux¬ 
iliary photon pair in a maximally entangled Bell state to 
achieve 75% success rate [27j. Although methods for pro¬ 
ducing single entangled pairs, such as biexciton decay in 
a quantum dots |42j . are becoming more efficient, sponta¬ 
neous parametric down-conversion (SPDC) remains the 
most widely available source of photon entanglement. 

In addition to studying accurate auxiliary states, we 
also analyze the case where the auxiliary state is pro¬ 
duced by a SPDC source. In this case, the |(^+) auxiliary 
Bell state can be approximated by pumping a nonlinear 
crystal to produce a multi-photon entangled state [43l - 
145] . The auxiliary state would then be given by 

OO 

n—0 

w{n, t) = y/n + 1 sech^r tanh^r 

\(i>t )34 = , + biblr 10)34 , 

n\y/n -I-1 

where the parameter r is the SPDC interaction param¬ 
eter and is related to mean photon number of the states 
being generated. The index n here represents the pair 
number so that terms |())^) are entangled photon states 
with 2n photons. 


III. METHODS 

A. Scheme quality: definitions and analysis 

An ideal destructive BSM scheme must unambigu¬ 
ously distinguish each of the Bell states without prior 
information on the state. The quality of a BSM scheme 
is described by the measurement success rate and fidelity. 

The maximum total success rate of a Bell-state ana¬ 
lyzer is defined as the unweighted average of the proba¬ 
bilities of unambiguously identifying each Bell state indi¬ 
vidually mms]. This standard definition for Bell-state 
analyzer success rate requires the input to be in the span 
of the four Bell states, which form a complete basis for 
the bi-photon Hilbert space. That is, by standard def¬ 
initions, a Bell-state analyzer performs a measurement 
on a bi-photon state. Therefore, this method gives an 
accurate representation of a Bell-state analyzer’s perfor¬ 
mance if the input is well-approximated by a bi-photon 
state, otherwise the input is not in the span of the Bell 
basis and the Bell-state analyzer would no longer be pro¬ 
jecting states onto just the Bell basis. We aim to pro¬ 
vide a comparison of the two Bell-state analyzers that is 


relatively independent of the exact input photon state; 
hence, input states containing more than two photons 
are not treated in this work. These cases complicate the 
standard definitions of a Bell-state analyzer, often re¬ 
quiring application-specific analysis to compute success 
rates. 

The maximum total success rate represents a 
quantum-limited upper bound on the efficiency of a Bell- 
state analyzer. The actual success rate of the measure¬ 
ment is less than the maximum because it depends on 
the choice of measurement outcomes used to identify 
the state and includes effects caused by imperfections 
in equipment. To illustrate this point, we consider the 
standard scheme. Applying the beam splitter transfor¬ 
mation given by Eq. © to each Bell state results in the 
output states 


\ip+) ^ (albl+albl'j |0) 





where was replaced with or b^ appropriately. For 
a perfect PNRD detector model, each output state is 
faithfully reported by the detector and so we are able 
to distinguish outputs of \ip^) from each other and from 
the outputs of 10^). However, both the \4>^} outputs are 
identical but for phase, thus the detector cannot distin¬ 
guish between them. Since only half of the Bell states are 
uniquely identified by measurement outcomes, we have 
a scheme maximum success rate of 50%. 

The BD model cannot resolve photon number and re¬ 
ports states aj^ |0) and |0) as a| |0) and |0) respec¬ 
tively. Fortunately these are still distinguishable from 
the outputs, again giving 50%. 

Our slow PNRD model reports states aJSJ |0) as aj |0) 
due to photons created by arriving during detector 
dead time. This change only affects the outputs of 
causing them to appear as one-photon events. Luckily, 
number resolution allows the detector to faithfully rep¬ 
resent the two-photon outputs of \(j)^) and thus can 
still be distinguished from !</)=*=) based on photon count. 
Since the outputs of !'(/'“) are also unaffected, this detec¬ 
tor model still allows for a 50% maximum success rate. 

The slow BD model, however, cannot faithfully re¬ 
solve two-photon states and hence cannot distinguish 
between output states of and |(^^), resulting in a 
detector-limited 25% maximum success rate. This 25%- 
efficient Bell measurement is often performed using time- 
bin-encoded photons if fast detectors are unavailable im- 

ssj. 

The enhanced scheme can be analyzed in a similar 
manner as the standard scheme but, for brevity, this will 
not be shown since the \ip^) Bell states have 80 possible 
outputs each, |(^+) has 42, and \(j)~) has 38. The unam¬ 
biguous output states of the enhanced scheme include 
four-mode, three-mode, and two-mode occupied states. 
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All single-mode states, where all four photons occupy the 
same mode, are ambiguous and must be rejected. Since 
there exists four-mode output states, where all photons 
occupy different modes, the success rate is nonzero when 
using non-number resolving detectors. The success rates 
for the enhanced scheme are given in section |IV A[ 

In addition to reduced success rates, including im¬ 
perfections in the analysis can lead to a reduction in 
the measurement discrimination fidelity of the scheme. 
To quantify this, we organize the measurements into 
two mutually exclusive categories: the positive (post- 
selected) measurements, and the negative (rejected) 
measurements. Within the positive measurements, a re¬ 
sult can either be true (successful) or false (unsuccess¬ 
ful). Hence, true positive measurements are those that 
collapse the state into a Bell state and correctly iden¬ 
tify which Bell state it collapsed into. False positive 
measurements are those that are post-selected but ei¬ 
ther do not collapse the state into a Bell state or falsely 
identify which Bell state it collapsed into. In this way, 
the true success rate of a BSM can be described by the 
true-positive rate, denoted here by pt ■ If we denote the 
false-positive rate by p/, the total probability of a posi¬ 
tive measurement is p* -|-p/ and the fidelity / of a BSM 
can then be computed by 


/ 


Pt 

Pt +Pf 


( 6 ) 


B. Dealing with large system sizes 


For a small apparatus like the standard scheme, mea¬ 
surement outcomes and probabilities can be computed 
by hand. However, when enhanced schemes with addi¬ 
tional auxiliary states are analyzed, the Hilbert space 
can become very large and numerical methods become 
necessary. The size of an n-photon Hilbert space with k 
states per photon is 


Nh 


k + n — 1 
n 


{k + n — 1)! 
(fc — l)!n! 


(7) 


For example, the size of the Hilbert space for the stan¬ 
dard scheme (n = 2 and fc = 4) is 10. However, the 
4-detector enhanced scheme (n = 4 and k = 8) has a 
Hilbert space size 330. Furthermore, if input and aux¬ 
iliary loss is accounted for, then the incident states can 
consist of zero, one, or two photons. In this case the size 
of the Hilbert space for the standard scheme increases to 
15 whereas for the enhanced scheme, it increases to 495. 

For number-resolving detector types, the number of 
measurement outcomes can be much larger than the 
Hilbert space size when including dark counts. For the 
enhanced scheme with PNRDs including detector inef¬ 
ficiencies, dark counts, and two auxiliary photons, the 
number of measurement outcomes is already 23392. The 
size of the Hilbert space exponentially increases when 
additional resources are added, with SPDC cases and 


some detector array numerical calculations in this study 
reaching Hilbert space sizes on the order of millions. For 
example, the largest auxiliary state we compute con¬ 
tains 20 photons, resulting in a Hilbert space size of 
Nh = 1560780 as determined using Eq. Q with k = 8 
and n = 22. 

Approaching these large computations using stan¬ 
dard positive-operator valued measure (POVM) meth¬ 
ods with the density operator formalism would require 
numerical arrays of sizes that can reach magnitudes 
of ~ 10^^ in some cases, which is impractically large 
for numerical methods. In addition, since the number 
of measurement outcomes can exceed the Hilbert space 
size, POVM operators can contain more than Nfj el¬ 
ements. For this reason, we approach the numerical 
problem by representing states as pure states rather than 
mixed states. This requires tracking the loss modes until 
measurement and treating dark counts as variations on 
measurement outcomes rather than modifications to the 
state itself. We expand our discussion on these methods 
in the following section. 


C. Loss, detector inefficiency, and dark counts 

To account for input photon loss prior to BSM, each 
Bell state was first passed through a beam splitter with 
weight y/ffi for input state transmission rate r]i. This 
transformation is given by 

+ \/l - ^7* > (8) 

where r is the photon creation operator for the loss 
mode. For the enhanced scheme, we applied an identical 
transformation on the auxiliary state using transmission 
rate parameter rja- This causes the output state prob¬ 
abilities to become dependent on iji and also ija for the 
enhanced scheme. 

To account for detector inefficiency, dark counts, and 
different detector models, we determined the condi¬ 
tional probabilities P(rni\ipj) that a measurement out¬ 
come rrii is triggered by a state \pj), where \ipj) form 
a basis for the scheme output Hilbert space. The to¬ 
tal probability that measurement outcome rrii occurs 
is then given by P{mi) = P(jni\ipj)P{ipj). This 

definition is consistent with a POVM {F^}, defined by 
Fi = j P{'mi\Lpj) \ipj) {pj\ so that the probability of 
measurement outcome rrii is P{mi) = Tr^Fip) where 
P = J2j,kPjk \Pj) {Pk\ and pjj = P{pj). Since by def¬ 
inition J2i = 1) we have J^i^t — 7- By de¬ 

scribing the measurement using conditional probabilities 
directly, we avoid the density operator formalism, which 
saves computational resources. 

Summing P(ini) over all post-selected measurements 
and averaging over each Bell state input gives the total 
probability of a positive measurement, pt+pf. The true 
success rate pt is determined in the same way with the 
exception that the sum is limited to unambiguous output 
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states only, as determined using = 1. The fidelity is 
then determined by the ratio of these sums. 

To calculate the conditional probabilities P{mi\(pj), 
each mode of each \ipj) of the output Hilbert space was 
first passed through a beam splitter with weight y/ffd for 
sub-unity detector efficiency ry^. This transformation is 
given by Eq. Q where rji is replaced by ry^;. By mod¬ 
eling detector inefficiency using a beam-splitter model, 
we can consider the parameter r]d as an effective effi¬ 
ciency/transmission parameter that implicitly accounts 
for any photon loss experienced between the scheme and 
the detectors, such as coupling loss. 

We then computed the probabilities associated with 
each outcome. At this point the loss modes were safely 
discarded. To include dark counts, we then allowed each 
mode to gain one photon with probability ^ or remain 
unchanged with probability 1 — Here we assume that 
^ <C 1 so that the probability of multiple dark counts 
occurring simultaneously in a single mode is negligible. 

If detectors other than PNRDs were being analyzed, 
each measurement outcome was modified according to 
detector model, either by eliminating photon count in¬ 
formation, eliminating the late mode if the early mode 
is detected, or both. This can be done only at the level 
of the probabilities, otherwise undesired interference can 
occur and lead to incorrect results. This method allowed 
us to simultaneously obtain probabilities P{mi\ipj) and 
identify the subset of relevant measurement outcomes. 


D. Detector array analysis 

We now proceed to describe the method we used to cal¬ 
culate analytic expressions for success rates using arrays 
of BDs. The BD-type models that we considered can be 
assumed to perform ideally if a single photon is incident 
when ignoring detector inefficiencies and dark counts. 
Therefore by splitting output modes using beam split¬ 
ters we can on average reduce the multi-photon detec¬ 
tion events that occur and recapture success rate lost by 
non-ideal detector types. The method we outline here to 
determine success rates is similar to the analysis method 
proposed by Kok and Braunstein m- 

We first give an example for the standard scheme and 
then extend this to include the enhanced scheme. To 
compute the success rate with increasing numbers of de¬ 
tectors, it is beneficial to first categorize output states 
based on the exponents of their photon creation oper¬ 
ators. For the standard scheme the two possible cate¬ 
gories are { 1 , 1 } and { 2 }, which for example would in¬ 
clude states such as a|S| and a|^, respectively. 

For the beam splitter detector array model, we always 
assume that one input mode of the beam splitter is the 
vacuum state with the other input being an output mode 
of the scheme. In this case, no interference can occur 
and so we can treat the array problem using classical 
statistics. 

Let P{1, 1} and P{2} denote the probabilities that an 


output of type (1,1} and {2} occur, respectively. These 
probabilities are computed assuming the input state is 
an equal distribution of each Bell state so that P rep¬ 
resents the average probability over the four Bell states. 
Similarly, let Pt{l, 1} and Pt{2} denote the probabilities 
that output states of those types also cause collapse into 
a Bell state, hence give rise to true-positive events. In 
this sense P{1, 1} -|- P{2} = 1 and Pt{l, 1} -I- Pt{2} = 5 . 
Although these sums always remain constant, in general 
the probabilities P and Pt can change with increasing 
detector number per array, denoted as N. In particular, 
since adding beam splitters can only ever convert type 
{ 2 } into type { 1 , 1 } states, probability must always flow 
from P{2} to P{1, 1} and from Pt{2} to Ht{l, 1}. 

To simplify notation, we introduce the probability vec¬ 
tors P = (P{1, 1 }, P{2})T and Pt = (Pt{l,l},Pt{2}r. 
For example, the standard scheme with BDs gives vec¬ 
tors P = ( 5 ,1)"’’ and Pt = ( 5 , 0 )'^ since all outputs of 
are of type { 1 , 1 } and unambiguous, whereas out¬ 
puts of 10^) are of type {2} and ambiguous (see Eq. (§)■ 

The linear transformation that governs the probability 
flow for the standard scheme after log 2 (N) beam splitters 
is given in the {{ 1 , 1 }, { 2 }} basis as 


.4(») = („ 1) 



(9) 


which can be obtained by analyzing the beam splitter 
transformation —)■ (c| +ic\)/\^ applied to each mode 
of archetype states representing {1,1} and {2}. Since 
BDs can only discriminate terms of type {1,1}, the prob¬ 
ability Pt{l,l} after applying A{N) is the true success 
rate. For example, for Pt = (|,0)"'’, we have that the 
true success ratept(iV) = {A{N)Pt)^i ij = ^ is constant 
and equal to the scheme maximum rate, as expected for 
the standard scheme with BDs. 

To extend this idea to the enhanced scheme, we consid¬ 
ered probabilities associated with five categories, namely 
{1,1,1,1}, {1,1,2}, {1, 3}, {2, 2}, and {4}. The trans¬ 
formation A{N) governing the probability flow in the 
{{1,1,1,1}, {1, 1 , 2 }, {1,3}, { 2 , 2 }, {4}} basis is given by 


(1 

N-1 

(Ar-2)(Af-l) 


(N-3)(N-2)(N-1) 

\ 

N 

JV2 

N2 

M3 



0 

1 

N 

3(Ar-i) 

N2 

2{N-1) 

N2 

6{N-2){N-1) 

M3 



0 

0 

1 

0 

4{N-1) 

N3 


.(10) 

0 

0 

0 

1 

Ar2 

3{N-1) 

N3 



VO 

0 

0 

0 

1 

ATS 




The true success rate is then given by (A(iV)Pt){i ^ ^ 
where Pt is obtained from the enhanced scheme outputs. 


If photon loss is taken into account, the number of 
categories for both schemes increases. For the standard 
scheme we must consider the additional two categories: 
{0} and {!}. Likewise, the enhanced scheme has seven 
additional categories: { 0 }, {!}, { 1 , 1 }, { 1 , 1 , 1 }, { 1 , 2 }, 
{2}, and {3}. 
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TABLE I. A comparison of the maximum success rates pt,max 
obtainable with unity fidelity for four different detector mod¬ 
els and two different schemes. These success rates are com¬ 
puted under the ideal case; with unity detector efficiency and 
photon transmission rjd = rji = rja = 1, and zero dark count 
probability ^ = 0. Note that the slow-prefixed detector mod¬ 
els are relevant to time-bin encoded applications only. 


Success Rate 

PNRD 

BD 

Slow PNRD 

Slow BD 

Standard Scheme 

1/2 

1/2 

1/2 

1/4 

Enhanced Scheme 

3/4 

3/16 

39/64 

1/16 


IV. RESULTS 

A. Scheme efficiency and fidelity 

In this section, we report success rates and fidelity for 
both schemes. For the enhanced scheme, we assume that 
the auxiliary state is produced by a single pair source 
that is subject to loss, representing pair generation inef¬ 
ficiency or propagation/coupling loss. 


1. Numerical results 


Firstly, without including loss, inefficiencies, or dark 
counts, we found that each detector type provides 
enough information to obtain a non-zero success rate 
with unity fidelity [TABLE |T]. For BD-type models, the 
optimal measurements to post select are unambiguous 
measurements that indicate each photon hit a different 
detector. The probability of this occurring was found 
to be less in the enhanced scheme than in the standard 
scheme. Thus, although the enhanced scheme outper¬ 
formed the standard scheme when PNRD-type models 
were used, it was inferior to the standard scheme for 
BD-type models. 

For realistic detector parameters [2H1 HH], the dark 
count probability can be between 10“"* and 10“®. For 
these probabilities, dark counts have negligible impact 
on success rates and fidelity for both Bell analyzers un¬ 
less detector efficiency is very small {r]d < 0.01). For the 
purposes of our numerics, we select a dark count proba¬ 
bility of ^ = 10"^ 

In repeater schemes. Bell-state analyzers are used to 
measure photons that must first travel over long dis¬ 
tances m This can cause photon losses around 99% 
before Bell-state analysis, effectively enhancing the in¬ 
fidelity caused by dark counts. When including input 
state loss, acceptable measurement hdelity can only be 
obtained when ^ ^ r]i,r]d. For this reason, we only ex¬ 
plore the regime where first-order dark counts dominate 
the false-positive rate. For the purposes of our numer¬ 
ics, we select an input state transmission of rji = 0.01, 
which is a realistic value in the context of quantum re¬ 
peaters. Since a successful Bell-state measurement only 
occurs when both input photons arrive at the analyzer, 








Detector efficiency, r]d 


Detector efficiency, rjj 


FIG. 2. (Color online) The effect of detector efficiency 
on the true success rate and fidelity for both the standard 
and enhanced Bell-state analyzers. Dashed and solid curves 
represent analytic approximations to the standard and en¬ 
hanced scheme solutions, respectively. Circle and diamond 
points represent numerically exact values for the standard 
and enhanced schemes, respectively. All plots are calculated 
with a dark count probability of ^ = 10“®, an auxiliary state 
transmission of r/a = 1.0, and an input state transmission of 
rji = 0.01. The success rates are plotted after dividing out the 
scheme-independent quantity rj^. The slow detector models 
are relevant to time-bin encoded applications only. 


the true success rate scales as pt oc r]f. This scaling is 
independent of the analyzer scheme and is divided out 
for a more direct scheme comparison [FIG.j^. 

Under dark count and loss conditions, our numerical 
results showed that the enhanced scheme with PNRDs 
provided higher success rates than the standard scheme 
when r]ii > 0.82, assuming rja = 1-0, or when rja > 0.82, 
assuming r]d = 1.0 [FIG. [^. For slow PNRDs, the en¬ 
hanced scheme can be superior when rid > 0.94, assum¬ 
ing rja = 1-0, or when r]a > 0.94, assuming pd = 1-0. 
For BD-types, the enhanced scheme showed no advan- 






















FIG. 3. (Color online) A closer view of the snccess rate inter¬ 
section points between the standard and enhanced schemes 
for PNRDs and slow PNRDs for different anxiliary state 
transmission values rja- Dashed and solid curves represent an¬ 
alytic approximations to the standard and enhanced scheme 
solutions, respectively. Circle and diamond points repre¬ 
sent numerically exact values for the standard and enhanced 
schemes, respectively. All numerical points are calculated 
with a dark count probability of ^ = 10 “® and an input trans¬ 
mission of rji = 0.01. The success rates are plotted after di¬ 
viding out the scheme-independent quantity 77 ^. The fidelity 
plots for different rja are not shown since, for rja S> they 
are nearly indistinguishable from those shown in FIG.[^ The 
slow PNRD detector model is relevant to time-bin encoded 
applications only. 


tage over the standard scheme for all efficiency values. In 
addition, the standard scheme was found to have better 
fidelity scaling with detector efficiency in all cases. 

It is important to note that all but the slow PNRD 
cases approach the scheme maximum rates given in 
TABLE H] as 77 ^ —>■ 1. We explain this discrepancy now. 
For the detector models labeled PNRD, BD, and slow 
BD, the ideal post-selected measurements do not become 
ambiguous when only loss and detector inefficiency are 
introduced. For these three detector types, it is possible 
to always identify and reject measurements that indi¬ 
cate fewer photons than expected. For PNRDs, this is 
possible due to ideal number-resolving capabilities. For 
BD-types, this is possible because the ideal post-selected 
measurements all already indicate photons arriving at 
different detectors. Therefore false-positive events can¬ 
not occur unless a false state has lost photons that are 
replaced by an equal number of dark counts during de¬ 
tection, which could then result in a positive measure¬ 
ment. Since dark counts are considered rare, the fidelity 
is already very good when detector efficiency is large for 
PNRD, BD, and slow BD models and so we need not 
reject additional measurements to improve fidelity. 

The slow PNRDs main advantage over the slow BD 
is its ability to distinguish states based on the photon 
count in the early mode. As shown in section [ill A| with 
the standard scheme, the slow PNRD allowed for a suc¬ 
cess rate of 50%, as opposed to the 25% maximum using 
a slow BD, by distinguishing from based only 
photon count information. As a consequence, when us¬ 
ing slow PNRDs, the fidelity is extremely susceptible to 
detector inefficiency and photon loss, even without the 
presence of dark counts. In the standard scheme, for 


example, detector inefficiency or photon loss can cause 
the two-photon states from \(j)^) to often appear as a 
single early photon, which are then falsely post-selected 
as indicating the Bell state. Thus, when using pho¬ 
ton count information from slow PNRDs, photon loss or 
detector inefficiency can cause false-positive events even 
in the absence of dark counts, if these specific types of 
measurements are post-selected. 

In order to maintain acceptable fidelity when using 
slow PNRDs, it is necessary to further reject measure¬ 
ments that become ambiguous when losses are intro¬ 
duced. This is accomplished by only accepting measure¬ 
ments that preserve the expected total photon count. 
For the standard scheme, this is equivalent to completely 
ignoring photon count information, effectively treating 
the slow PNRD as a slow BD, giving a maximum suc¬ 
cess rate of 1/4. For the enhanced scheme, this resulted 
in an improved success rate over the corresponding slow 
BD case with a maximum success rate of 9/32. These 
two values are the success rates for 77 ^ = 1.0 illustrated 
in the slow PNRD case in FIG. [21 


2. Analytic approximations 

Once the maximum success rates pt,max are deter¬ 
mined for a given detector model and scheme, under¬ 
standing the effect of losses and inefficiencies on the suc¬ 
cess rate is uncomplicated. The vast majority of true¬ 
positive measurements occur when detectors operate cor¬ 
rectly and no photons are lost. For the standard scheme, 
the two possible loss channels are input loss and detector 
inefficiency, represented by parameters rji and r]d, respec¬ 
tively. Since the standard scheme utilizes two incident 
photons, two detectors, and has four modes that expe¬ 
rience dark counts, the total success rate is given by 
Pt = VdPii^ - OV.max +P{, where pt.max is 1/2 in the 
PNRD and BD cases and 1/4 for the slow-type models. 
The small correction is a positive contribution to the 
success rate caused by true-positive cases resulting from 
inefficiencies and dark counts effectively canceling each 
other, leaving the measurement conclusion unaffected. 
For dark count probability ^ 1, we have (1 — ~ 1 

and ~ 0 leaving pt ~ rjjrj'^Pt.ma^- 

The enhanced scheme can be analyzed in a similar 
manner. In addition to the two loss channels listed 
above, the enhanced scheme can suffer from auxiliary 
state photon loss, parameterized by rja- When the en¬ 
hanced scheme has no auxiliary input, the scheme resorts 
to behaving like the standard scheme, with a maximum 
success rate of 50% m- Hence, when 77 ^ 7 ^ 1, some 
success rate could be recovered by post-selecting some 
measurements where both auxiliary photons are lost, in 
addition to post-selecting measurements where the aux¬ 
iliary state is unaffected. Unfortunately, if this is done, 
false-positive measurements can arise from input state 
photon loss directly, similar to the slow PNRD case dis¬ 
cussed previously. With large input loss, the reduction in 
fidelity caused by including measurements resulting from 
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a vacuum auxiliary state will not justify the gain in suc¬ 
cess rate, so we only consider post-selecting terms. 
Since the enhanced scheme utilizes four incident pho¬ 
tons, four detectors, and has eight modes that can ex¬ 
perience a dark count, the total success rate is given by 
Pt = wherept,max is 3/4, 3/16, 

9/32, and 1/16 for the PNRD, BD, slow PNRD, and 
slow BD cases, respectively. Again, for ^ <C 1 we have 
(1-^)® ~ 1 andp 4 ~ 0 so that pt ~ These 

zero-order terms in ^ well-approximate the numerical re¬ 
sults in FIG. [Hand FIG. [H 

With these results, it is straightforward to pre¬ 
dict the bounds where the enhanced scheme is supe¬ 
rior to the standard scheme. For P NRD s we have 
> v'i'nd/'^ implying r^dPa > \/^ ^ 0.82. For 
slow PNRDs we have 977 / 77 ^ 77^/32 > 77 / 77^/4 implying 
rjdila > 2-\/2/3 ~ 0.94. Both of these expressions agree 
with our numerical results. These results imply that 
the enhanced scheme is superior to the standard scheme 
in terms of success rate if elficient detectors with at 
least some number resolving capabilities are utilized and 
bright pair sources are available to produce the auxiliary 
state. 

Understanding the effect of loss and dark counts on 
fidelity can be done in a similar manner. Since only 
measurements that preserve the expected photon count 
are post-selected, all false-positive measurements must 
arise from at least one dark count event. In this work we 
assume ^ <C 77 ^ <C I so that only terms first-order in dark 
count probability dominate the false-positive rate. This 
assumption is justified in any regime where the fidelity 
is acceptable. It also provides an upper bound on the 
fidelity even in the regime ^^ 77 ^. In addition, if we 
further assume that rji ^ ild,Va, the dominant term will 
also correspond to the case where a single input photon is 
lost before arriving at the analyzer, and no other photons 
are lost. This assumption is reasonable for a repeater 
scheme, where input state transmission is poor. 

Knowing that, in the regime ^ ^ rji ^ %,? 7 a, the 
false-positive rate is dominated by one-photon input loss 
and single dark count events, the false-positive rates are 
of the form pf ~ Cir]drii{l — for the stan¬ 
dard scheme and pf ~ ^ 277 / 77 ^ 77^(1 — for the 

enhanced scheme. The constant prefactors, Ci and C 2 , 
depend on detector model and can be solved by directly 
counting the probability that single-photon inputs are 
post-selected, given the output states and post-selected 
measurements of each scheme. For example, for the stan¬ 
dard scheme with PNRD detectors, the post-selected 


measurements indicate either a\b\ | 0 ), | 0 ), a\b[ | 0 ), 


or a\h\ |0) (see section III A). The probability of obtain¬ 


ing a single photon output is 2777(1 — 777 ) for each mode. 
This single photon can only become post-selected if a 
dark count occurs in two of three possible modes. For 
example, a\ | 0 ) would require a dark count in either b\ 
or Sj to give rise to a false positive measurement. Hence, 
in this case, the prefactor is Ci = 4. Likewise, for the 


slow BD case, only measurements indicating ajSJ |0) or 


a\b\ |0) are post-selected. Here, a single photon can only 
become post-selected if a dark count occurs in one of 
three possible modes. Hence the prefactor is Ci = 2, 
where the 2 arises from the probability of obtaining a 
single photon output. 

The prefactors Ci for the standard scheme were found 
to be 4 for PNRD and BD models and 2 for the slow de¬ 
tectors. The enhanced scheme prefactors C 2 were found 
to be 10, 5/2, 8/3, and 3/4 for PNRD, BD, slow PNRD, 
and slow BD models, respectively. These false positive 
rates were used along with the analytic success rate ap¬ 
proximations to determine the fidelity. The resulting 
fidelity well-approximates the numerically exact results 
in FIG. [Hand FIG.! 

An interesting consequence of the result that both 
Pt oc 77 / and p/ oc 77 / is that the fidelity, pt/{pt +Pf), 
is independent of the auxiliary state photon loss rate. 
In addition, pf ex. ^ for ^ <C 777 ^ 1 indicates that any 
difference in fidelity between the standard and enhanced 
schemes can be made negligible with a small enough de¬ 
tector noise. Decreasing the dark count probability be¬ 
low 10 “® would show significant improvements to fidelity 
with inefficient detectors for both schemes. 


B. Detector arrays 

In this section, we report success rates and fidelity 
for both schemes assuming the detectors are replaced 
by arrays of BD or slow BD detectors. For the en¬ 
hanced scheme, we still assume that the auxiliary state 
is produced by a single pair source that is subject to 
loss, representing pair generation inefficiency or propa¬ 
gation/coupling loss. 


1. Numerical results 

Firstly, with high detector efficiency and a very bright 
auxiliary pair source, we found that arrays of BD or slow 
BD types can allow the enhanced scheme success rate to 
exceed that of the standard scheme [FIG.|^. Note that 
the success rate for the standard scheme using BDs is al¬ 
ready optimal without using detector arrays and adding 
arrays of BDs to the standard scheme would only de¬ 
crease the fidelity. For this reason, we compare the en¬ 
hanced scheme BD case to the optimal iV = 1 values for 
the standard scheme BD case. In this case, the enhanced 
scheme becomes superior to the standard scheme when 
four detectors are used per output (16 total). 

For slow BDs, using detector arrays improves the suc¬ 
cess rate for both the standard and enhanced schemes. 
However, due to a larger maximum success rate, the en¬ 
hanced scheme can surpass the standard scheme with ar¬ 
rays of at least five detectors per output (20 total). Note 
that a uniform array of size other than a power of two 
cannot be implemented using just 50/50 beam splitters 
but could be implemented using multi-port beam split¬ 
ters or uneven (non 50/50) beam splitters. This model 
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also gives upper bound rates for other spatial multiplex¬ 
ing methods such as beam-shaping optics [40]. 

Adding arrays of non-number-resolving detectors will 
only ever recover success rate up to the success rate 
achievable with PNRDs. For the enhanced scheme to be¬ 
come superior to the standard scheme using either BDs 
or slow BDs with increasing array size, we still require 
VaVd > 0.82 (see section IVA). Decreasing the detector 
efficiency or the auxiliary state transmission decreases 
the success rate, as expected, and moves the intersec¬ 
tion point between the standard and enhanced scheme 
success rate to larger array sizes. 

When using large arrays, the detector dark counts 
have a large impact on the fidelity when input state pho¬ 
ton loss is high, even if detector efficiency is high [FIG. 
1^. In the array size regime where the enhanced scheme 
overtakes the standard scheme in success rate {N > 4), 
a dark count probability of ^ < 10“® is required for ac¬ 
ceptable hdelity. Numerically computing the hdelity and 
success rate becomes computationally difficult for array 
sizes A^ > 8 for the enhanced scheme. This is partially 
due to the large Hilbert space size, which is 814385 for 
N = 8 when including photon loss, but also due to 
the very large number of measurement outcomes, which 
could be up to X]i=o (T) “ 8303633 in the BD case for 
4 photons and one dark count in 64 modes. Monte Carlo 
methods might prove a better method to compute nu¬ 
merical values for larger array sizes m- In the following 
section, we determine analytic approximations that give 
insight into the impact of imperfections on the success 
rate and fidelity for larger array sizes. 


2. Analytic approximations 




FIG. 4. (Color online) Increase in true success rates using ar¬ 
rays of imperfect detectors to approximate photon number re¬ 
solving detectors. Dashed and solid curves represent analytic 
approximations to the standard and enhanced scheme solu¬ 
tions, respectively. Circle and diamond points represent nu¬ 
merically exact values computed for computationally reason¬ 
able N. The fidelity curves for the enhanced scheme (solid) 
are lower bound approximations obtained from upper bound 
false-positive rates. The fidelity curves and data points were 
computed using input state transmission pi = 0.01, detector 
efficiency rjd = 1.0, two different dark count probabilities, 
^ = 10“® or ^ = 10“®, and an auxiliary state transmission 
r]a = 1.0 when relevant. The success rates are plotted af¬ 
ter dividing out the scheme-independent quantity rji. The 
dashed lines for the BD case show the N = 1 values for 
the standard scheme since arrays of BDs cannot improve the 
standard scheme success rate. The slow BD detector model 
is relevant to time-bin encoded applications only. 


To make analytic predictions for the success rate and 
the fidelity, we use the methods described in section 
[rag Recall that, as an example, in section mE we 
showed that the standard scheme with BDs is associated 
with probability vectors P = (^, 1)"^ and Pt = (|,0)'^ 
in the {{1,1}, {2}} basis. This resulted in the constant 
success rate pt{N) = 1/2. Now, we use the same reason¬ 
ing to solve the standard scheme with slow BDs as well 
as the enhanced scheme with BDs and slow BDs. 

To determine the success rate expression for the stan¬ 
dard scheme with slow BDs, we first found the effective 
probability vectors for the slow BD model. These vectors 
are different from the regular BD case and, at the level of 
the probabilities, can be computed from the scheme out¬ 
puts by combining early and late time-bin modes of the 
same spatial mode if the early mode is occupied. Notice 
that in this case we do not wish to eliminate the pho¬ 
tons in the late mode since for this method to work cor¬ 
rectly the total photon count must be preserved. That 
is, we need only capture the behavior of the ‘slow’ flaw 
since the non-number-resolving flaw of the slow BD is 
automatically captured by only allowing probability of 
{1,1} types to contribute to success rate. Performing 


this modification to the output states of the standard 
scheme gives the vectors P = (q, f)"'" and Pt = (j, j)"'"- 
By applying the array transformation A{N) to Pt we get 
the success rate 

PtiN) = (A(iV)P,){, = 1-^- 

The enhanced scheme output gives probability vectors 
describing the initial distribution over the categories in 
the basis {{1,1,1,1}, {1,1,2}, {1,3}, {2,2}, {4}} as 

pj^ 1 1 1 AY P =(^1 A A 1 oY (12) 

I 4’16’16’32’32 / ‘ V 16’16’16’16’ / 


for the BD model. Applying the transformation A{N) 
for the enhanced scheme gives success rate 


Pt{N) - (A(iV)A){ia,i.i} - 4 “ + PPPp ■ (13) 

Likewise, for the slow BD model we have 


P = 


/I 1 1 7 


V8’ 4’ 4’ 32’ 

32 / 


Pt = 


I I I I I 
16’ 4’ 4’ 8’ l6 


(14) 
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which results in 


Pt{N) 


3 13 , 21 3 

4 ~ ^ 161V2 “ ■ 


(15) 


Taking N ^ oo we see that pt {N) approaches the scheme 
maximum rates in all cases, as expected. 

To include dark counts, detector inefficiency, and pho¬ 
ton loss, we again assumed that the vast majority of 
true-positive events occur when detectors operate cor¬ 
rectly and all photons arrive at the analyzer. This gives 
Pt ~ ~ for the standard scheme and 

Pt ~ 1d''laVi{^ ~ for the enhanced scheme. 

The term AN and 8N are the total number of modes 
after applying an A^-sized beam splitter model for the 
standard and enhanced schemes, respectively. 

To estimate the false-positive rate while including im¬ 
perfections, we use the same arguments as in the N = 1 
case in section |1V A| However, for detector arrays, the 
prefactors Ci and C 2 become dependent on TV in a non¬ 
trivial way. For the standard scheme, the prefactor is un¬ 
complicated because most false-positive measurements 
arise from states of type {1}, and P{1} = 2r]i{l — rji) 
is independent of N. For slow BDs, by inspection, the 
number of false-positive combinations was determined to 
be (27V—1), and so C'i(7V) = 2(27V—1). Hence the false¬ 
positive rate for the standard scheme with arrays of TV 
slow BDs is approximately 


Pf ~ 2(27V - l)rjdV^{l - (16) 


a very slight decreasing trend with increasing TV, im¬ 
plying 2/3 can provide a more accurate upper bound 
false-positive rate. Thus the false-positive rate for the 
enhanced scheme with arrays of TV BDs is approximated 
by 


Pf 3 ) 773^(1 - 1,1}, (19) 

where P{1,1,1} is given by Eq. Note that when 

TV = 1, we recover the result from section |IV A| for BDs. 
That is, (2/3)(8TV-3)(2-2/Af-h3/4TV2) = 5/2 for = 1. 
Since the 2/3 factor is an upper bound, the fidelity de¬ 
rived from this false-positive approximation will slightly 
underestimate the numerical values for fidelity at higher 
N [FIG.|^. 

For slow BD detectors, {l,l,l}-type states have 
8TV — 6 modes that could suffer a dark count to become a 
{1,1,1,1} type. This is because a dark count must occur 
in a different spatial mode than those already occupied 
to change the total measured photon count, as a conse¬ 
quence of the slow recovery property. This gives an up¬ 
per bound Pf < (8TV — 1,1}. For 

the TV = 1 slow BD case, the fraction of rejected states 
can be directly counted to be 3/4, again by inspection 
of the scheme output states. This fraction was again 
found to have a decreasing trend, slightly more signifi¬ 
cant than in the BD case. Thus the false-positive rate 
for the enhanced scheme with arrays of TV slow BDs is 
approximately 


The enhanced scheme false-positive rates are more dif- 
hcult to determine because most false-positive measure¬ 
ments arise from states of type {1,1,1}, and P{1,1,1} 
increases with N due to contribution from states of type 
{1,2} and {3} splitting. Using A{N) with input state 
and auxiliary state loss considerations, the TV-dependent 
P{1,1,1} is 

Pjl, 1,1} ~ 772 , 7,(1 - , 7 ,) (^2 - ^ (17) 


Pf ^(8TV - 6 ) 773^(1 - 1,1}, (20) 

where P{1,1,1} is given by Eq. ( [I^ . With TV = 1, 
we again recover the result from section |IV A| for slow 
BDs. That is, (3/4)(8AT - 6)(2 - 5/2N +1/N^) = 3/4 
for TV = 1. The fidelity derived from this false-positive 
approximation will underestimate the numerical values 
for fidelity at higher TV more so than the corresponding 
BD fidelity [FIG.|4]. 


for the BD model, assuming rji ^ rja- Likewise, for the 
slow BD model we have 


C. Alternative auxiliary source 


P{l,l,l}c.77a^.(l-77,)(^2- A + (18) 

An approximate upper bound for the false-positive 
rate is given by the probability that a state of type 
{1,1,1} is measured to be of type (1,1,1,1}. Eor BD 
detectors, {l,l,l}-type states have 3 of the 8TV modes 
occupied and so 8TV — 3 modes could suffer a dark count 
to become a {1,1,1,1} type. This gives an upper bound 
Pf < (8TV — 3)77^^(1 — ^)®^“^P{1,1,1}. This is an up¬ 
per bound because a fraction of the (1,1,1,1} measure¬ 
ments are actually rejected during post-selection. Eor 
the TV = 1 BD case, the fraction of rejected states 
can be directly counted to be 2/3 by inspection of the 
scheme output states. This fraction was found to have 


We now briefly discuss the case where the enhanced 
scheme auxiliary state is generated by a SPDG source. 
In this section we still consider the input state |'I')j ^2 fo 
be a bi-photon state so that the state bei ng me asured is 
in the span of the Bell basis (see section HI A[ ). In this 
case, |4 >)j^2 \'Pt )34 has constant photon number equal to 
2n + 2. This implies that states with different n are or¬ 
thogonal and will not interfere when passing through the 
analyzer. Hence it suffices to analyze the performance 
of the scheme using auxiliary states \'Pt )34 separately 
and average the resulting rates using the weight func¬ 
tion w‘^{n,T) (see Eq. ©)■ 

If the input is a superposition of states with different 
photon numbers (such as produced by another SPDG 
source), then this analysis method fails. However, the 
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Auxiliary pair number, n 


PDC interaction parameter, T 


FIG. 5. (left) Success rates and fidelity for the enhanced 
scheme with different multi-photon auxiliary states from a 
SPDC source, (right) True success rates and hdelity, assum¬ 
ing a SPDC-based auxiliary source, after a weighted average 
of the true-positive and false-positive rates using weighting 
function w^{n, r). Curves on the right correspond to PNRD 
(solid), slow PNRD (dashed), BD (dotted), and slow BD 
(dot-dashed) detector types. All plots were computed in the 
ideal case; with unity detector efficiency and photon trans¬ 
mission rjd = rji = rja = I, and zero dark count probability 
^ = 0. The slow-prefixed detector models are relevant to 
time-bin encoded applications only. 


input state need not be a true bi-photon state to test 
the Bell-state analyzer. Two high-quality single photons 
could be heralded from weakly-pumped SPDC sources 
m and then projected onto a Bell state by the measure¬ 
ment. The efficiency and fidelity of the measurement can 
be inferred using a decoy state protocol, as if the input 
was a true bi-photon state [11 fi) . 

To analyze the enhanced scheme using a SPDC auxil¬ 
iary source, we begin by determining the maximum suc¬ 
cess rate for each detector type, in the absence of dark 
counts, detector inefficiencies, and loss rates. To do this, 
we consider each |(/)+) for 0 < n < 10 as an auxiliary 
state separately. For each n, the true-positive measure¬ 
ments were determined and then combined to form four 
mutually exclusive sets of measurements that are used 
to identify each of the four Bell states. These combined 
sets of post-selected measurements were then used to 
determine the true-positive and false-positive rates for 
each auxiliary state \(j)^) along with the associated fi¬ 
delity. This was done for each detector type, the results 
of which are summarized in FIC. [5] 

Interestingly, the maximum success rate for the PNRD 
case remains above 50% regardless of the auxiliary pair 
number. This implies that the enhanced scheme with 
a SPDC auxiliary source still has an upper bound rate 
greater than that of the ideal standard scheme. The 


maximum success rate for the enhanced scheme using 
a SPDC source is 58%, which occurs with an interac¬ 
tion parameter value of r = 0.67, corresponding to a 
mean photon pair number n = 1 for the auxiliary source. 
Values for mean photon pair numbers when generating 
Bell pairs using SPDC are generally chosen to be small 
(- 0.1) a, since the multi-photon states are undesired. 
However, even though the enhanced scheme operates ide¬ 
ally with a single Bell pair, FIG. [^suggests that the op¬ 
timal mean photon pair number is much larger than 0.1 
since it is possible to utilize both the empty pulses and 
the n > 1 entangled states generated by SPDC to per¬ 
form a successful BSM. If the distribution of pair num¬ 
bers can be focused nearer to n = 1, it may be possi¬ 
ble to raise the success rate higher than 58%. Perhaps 
this could be achieved by antibunching emission of pho¬ 
ton pairs from the SPDC source using a quantum Zeno 
blockade (52] , although this could never increase the rate 
above the 75% maximum. 


Since PNRDs can indicate photon number, each mea¬ 
surement can be attributed directly to one of the terms 
in the auxiliary state. For this reason combining unam¬ 
biguous measurements from each n results in a set of 
measurements that give perfect hdelity for any interac¬ 
tion parameter. Similarly, the partial number-resolving 
capability of the slow PNRD model allows a set of unam¬ 
biguous measurements to exist, which allows measure¬ 
ments with perfect Hdelity, albeit at an inferior rate to 
the standard scheme. On the other hand, BD-types can¬ 
not indicate photon count and so measurements from one 
pair number can be confused with those from another, 
leading to unavoidable false positive events. 


Computing success rates and Hdelity for large auxiliary 
states is computationally difficult even without including 
losses, inefficiencies, and detector dark counts. For n = 
10 , the computation involves 22 photons, leading to an 
output Hilbert space of size Nh = 1560780 (see section 


HI B). Including all photon losses increases the size to 


5852925. Since the number of measurement outcomes 
can be even larger than the Hilbert space size when using 
PNRDs with detector dark counts, we do not attempt 
numerical computations accounting for these additional 
imperfections. Instead, we make some observations for 
small n to make a deHnitive comparison between the two 
schemes when using an SPDC auxiliary source. 


For the enhanced scheme to be superior to the stan¬ 
dard scheme when using PNRDs with an SPDC auxil¬ 
iary source, it is necessary to accept many measurements 
that indicate different n. As a consequence, the proba¬ 
bility of losing two photons can directly contribute to 
inHdelity in the absence of dark counts, similar to the 
slow PNRD case discussed in section HV A1 Consider the 
incident states | 00)^2 1 ^ 2)34 Due to 

the symmetry of the apparatus, these two states produce 
many identical outputs. Unsurprisingly, we found that 
| 00 )i 2 |^J )34 can produce all of the measurement results 
that are post-selected to indicate a |<()''') 12 I<(' 1 ') input. 
This implies that pf oc {1 — assuming rji ^ Pditla 




































13 


and ^ <C 1. Thus for small rji, the fidelity is likely 
also small, unlike previous cases where the false-positive 
measurements were dominated by single detector dark 
count cases. For this reason, SPDC sources are not a 
suitable replacement for the auxiliary source in the en¬ 
hanced scheme when rji 1. This alternative auxil¬ 
iary source could still be used in a local application with 
very efficient number resolving detectors and low pho¬ 
ton loss rates. In this case, the results shown in FIG. 

become relevant. Regardless, an increase of 0.08 in 
maximum success rate over the standard scheme is likely 
not enough to justify the required increase in apparatus 
quality and complexity. 

V. CONCLUSIONS 

In this paper, we determined realistic success rates and 
measurement discrimination fidelities expected when im¬ 
plementing an enhanced linear-optic Bell measurement 
scheme using non-ideal detectors and incident states. 
This enhanced scheme was compared to the standard 
scheme for implementing a Bell-state analyzer. We de¬ 
scribed the numeric methods that we used to quantify 
the quality of the Bell-state analyzers in light of the chal¬ 
lenges caused by large Hilbert space sizes when detector 
arrays and complicated auxiliary states are considered. 
In addition, we used analytic methods to determine ap¬ 
proximate expressions that give insight into the effect of 
different imperfections on the scheme quality. 

When detector inefficiencies, dark counts, input pho¬ 
ton losses, and auxiliary photon losses were accounted 
for, the enhanced scheme using PNRDs was found to 
be superior to the standard scheme when detector and 
auxiliary source efficiency was high, 77^770 > 0.82. When 
using single non-PNRDs, the enhanced scheme showed 
no advantage over the standard scheme. However, the 
enhanced scheme’s 75% maximum rate allowed it to sur¬ 
pass the standard method when used with efficient non- 
PNRDs in array configurations. 

Both schemes that were analyzed suffered decreases in 
success rate when used with time-bin qubits and detec¬ 
tors with long dead time. In this case, arrays of highly 


efficient detectors again improved success rates arbitrar¬ 
ily close to scheme maximums in the large array limit. 
However, detector arrays are significantly susceptible to 
detector dark counts since the chance of misfire during 
a detection window is much higher with more active de¬ 
tectors. This effect is compounded by input loss, causing 
a reducing fidelity with increasing array sizes. If detec¬ 
tors composing the arrays have dark count probabilities 
smaller than 10 “®, the measurement fidelity can be made 
acceptable even with large input loss. 

We analyzed whether a parametric down-conversion 
source can be used to produce the auxiliary photon state 
required for the enhanced scheme. With this alternative 
auxiliary source and number-resolving detectors, the en¬ 
hanced scheme had a maximum success rate of 58%, still 
greater than that of the standard scheme. However, the 
measurement discrimination fidelity is very low when in¬ 
put loss is large. Hence, using a SPDC source with the 
enhanced scheme would only be advantageous for local 
applications with small loss rates and efficient number¬ 
resolving detectors. 

Although the regime where this enhanced scheme sur¬ 
passes the standard scheme is narrow, the advent of ac¬ 
curate, bright pair sources and improved detectors might 
allow an auxiliary photon assisted Bell-state analyzer 
to become viable in the near future. Furthermore, the 
methods and results detailed in this paper may prove 
useful for the analysis of the performance of other linear 
optical quantum information processing protocols under 
realistic conditions. 
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